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The possibilities of unlike particle correlations for a study of the space-time asym- 
metries in particle production, including the sequence of particle emission, are 
demonstrated. 



1 Introduction 

The correlations of particles at small relative velocities are widely used to study 
space-time characteristics of the production processes. Particularly, for non- 
interacting identical particles, like photons, this technique is called intensity or 
particle interferometry. In this case the correlations appear solely due to the 
effect of quantum statistics (QS).^'^ Similar effect was first used in astronomy 
to measure the angular radii of stars by studying the dependence of the two- 
photon coincidence rate on the distance between the detectors (HBT effect 
■^). In particle physics the QS interference was first observed as an enhanced 
production of the pairs of identical pions with small opening angles (GGLP 
effect ^). Later on, similar to astronomy, Kopylov and Podgoretsky^ suggested 
to study the interference effect in terms of the correlation function. ° 

The effect of QS is usually considered in the limit of a low phase-space 
density such that the possible multi-particle effects can be neglected. This 
approximation seems to be justified by present experimental data which does 
not point to any spectacular multi-boson effects neither in single-boson spectra 
nor in two-boson correlations. These effects can be however of some importance 
for realistic simulations of heavy ion collisions.'* They can also clearly manifest 
themselves in some rare events (e.g., those with large pion multiplicities) or 
in the eventually overpopulated regions of momentum space (see, e.g.f'^ and 
references therein).'' 



"Note that though both the KP and HBT methods are based on the QS interference, they 
represent just orthogonal measurements.^ The former, being the momentum-energy mea- 
surement, yields the space-time picture of the source, while the latter does the opposite. In 
particular, the HBT method provides the information about the angular size of a star but, 
of course, - no information about the star radius or its lifetime. 

''It was shown,^'*' using the analytically solvable model with the Gaussian distribution of the 
emission points characterized by the dispersions and in the momentum and ordinary 
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The particle correlations are also influenced by the effect of particle inter- 
action in the final state (FSI)7'* Thus the effect of the Coulomb interaction 
dominates the correlations of charged particles at very small relative momenta 
(of the order of the inverse Bohr radius of the two-particle system), respec- 
tively suppressing or enhancing the production of particles with like or unlike 
charges. Regarding the effect of the strong FSI, it is quite small for pions, 
while for nucleons it is often a dominant one due to the very large magnitude 
of the s-wave singlet scattering length of about 20 fm. 

Though the FSI effect complicates the correlation analysis, it is an impor- 
tant source of information allowing one to measure the space-time characteris- 
tics of the production process even with the help of non-identical particles?'^ 
Moreover, the unlike particle correlations, contrary to those of identical parti- 
cles, are sensitive to the relative space-time asymmetries in their production, 
e.g. - to the relative time delays, thus giving an important complementary 
information not accessible to the standard interferometry measurements.^" In 
the following we will briefly formulate the theory of these correlations and 
demonstrate the possibilities of the corresponding correlation technique. 

2 Formalism 

As usual, we will assume sufSciently small phase-space density of the produced 
multi-particle system, such that the correlation of two particles ciniitted with a 
small relative velocity in nearby space-time points is influenced by the effects 
of their mutual QS and FSI only.'^ We define the ideal two-particle correlation 
function R{pi,p2) as a ratio of the differential two-particle production cross 
section to the reference one which would be observed in the absence of the 
effects of QS and FSI. In heavy ion or high energy hadronic collisions we can 
neglect kinematic constraints and most of the dynamical correlations and con- 
struct the reference distribution by mixing the particles from different events. 

Assuming the momentum dependence of the one-particle emission prob- 
abilities inessential when varying the particle 4-momenta pi and p2 by the 
amount characteristic for the correlation due to QS and FSI {smoothness as- 
sumption), i.e. assuming that the components of the mean space-time distance 
between particle emitters are much larger than those of the space-time extent of 
the emitters, we get the well-known result of Kopylov and Podgoretsky for iden- 
tical particles, modified by the substitution of the plane wave e'^^^i+'P^^^ by 

space respectively, that a pion condensate of a small momentum dispersion of^ A/2ro <C 
gradually develops with the increasing phase-spa<;e density. 

■^This assumption may be not valid in the case of low energy heavy ion reactions when the 
particles are produced in a strong Coulomb field of residual nuclei. To deal with this field a 
quantum adiabatic approach can be used.^" 
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the nonsymmetrized Bethe-Salpeter amplitudes in the continuous spectrum of 
the two-particle states i]jp[p}{xi, X2), where Xi = {ti, r^} are the 4-coordinates 
of the emission points of the two particles and S is their total spin? At equal 
emission times in the two-particle c.m.s. {t* = t\ — t\ =0) this amplitude 
coincides (up to an unimportant phase factor due to the c.m.s. motion) with a 
stationary solution of the scattering problem '^.J^. (r*), where k* = p| = — 
and r* = — (the minus sign of the vector k* corresponds to the reverse 
in time direction of the emission process). The Bethe-Salpeter amplitude can 
be usually substituted by this solution {equal time approximation)."* Then, for 
nonidentical particles, 

R{Pi,V2) = Y.Ps{\&X^*)?)s. (1) 
s 

Here the averaging is done over the emission points of the two particles in a 
state with total spin S populated with the probability ps, '^s PS ^ ^-"^ ^'^^ 
identical particles, the amplitude in Eq. (1) should be properly symmetrized: 

V-flt^Cr*) - [V!Lt^(r*) + (-l)^Vf^^(r*)]/^/2. (2) 

3 Measuring the relative space-time asymmetries 

The correlation function of two nonidentical particles, compared with the iden- 
tical ones, contains a principally new piece of information on the relative space- 
time asymmetries in particle emission.'^'^ This is clearly seen in the case of 
neutral particles when the two-particle amplitude tp_l^. {r*) takes on the form 

V!Lt)(r*)=e-'''*-*+<^f.(r*), (3) 

where the scattered wave (/>f. (r*), in the considered region of small relative 
momenta, is independent of the directions of the vectors k* and r* . Inserting 
Eq. (3) into the formula (1) for the correlation function, we can see that the 
latter is sensitive to the relative space-time asymmetry due to the odd term 

''The equal time approximation is valid on condition* \t*\ <C m,2,ir*^ for sign(t*) = ±1 
respectively. This condition is usually satisfied for heavy particles like kaons or nucleons. But 
even for pions, the t* = approximation merely leads to a slight overestimation (typically 
< 5%) of the strong FSI effect and, it doesn't influence the leading zero-distance {r* <S \a\) 
effect of the Coulomb FSI. 

'^For unpolarized particles with spins si and S2 the probability pg = (25-|-l)/[(2si + I)(2s2 + 
1)]. Generally, the correlation function is sensitive to particle polarization. For example, if 
two spin-1/2 particles arc emitted with polarizations Pi and P2 then* po = (1 — Pi • P2)/4 
and pi = (3 + Pi • P2)/4. 
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^ sink*r*. Particularly, it allows for a measurement of the mean relative 

delays (t) = {ti — 12) in particle emission. To sec this, let us make the Lorcntz 
transformation from the rest frame of the source to the c.m.s. of the two 
particles: = ^{tl — vt), = Tt- Considering, for simplicity, the behavior 
of the vector r* in the limit \vt\ 3> r, we see that this vector is only slightly 
affected by averaging over the spatial distance r <C \vt\ of the emission points 
in the rest frame of the source so that r* —'yvt. Therefore, the vector r* is 
nearly parallel or antiparallel to the velocity vector v of the pair, depending 
on the sign of the time difference t = At = ti — 12- The sensitivity to this sign 
is transferred to the correlation function through the odd in k*r* w —'yk*vt 
term provided the sign of the scalar product k*v is fixed. 

For charged particles there arise additional odd terms due to the conflu- 
ent hypergeometrical function F{a, 1, 2;) = 1 + az/V? + a{a + l)z'^/2\^ + ■ . ., 
modifying the plane wave in Eq. (3): 

&J{rn = e''./Aj^[e-'''''''F{-iv,l,ip) + cl>f,,{r*)] , (4) 

where p = k*r* + k*r*, ij = {k*a)~^, a = {pziZ2e'^)~^ is the Bohr radius of 
the two-particle system taking into account the sign of the interaction (zie 
are the particle electric charges, fi is their reduced mass), S = argr(l + jt/) 
is the Coulomb s-wave shift and Ac{t]) = 2ttt} / [exp{2nr]) — 1] is the Coulomb 
penetration factor.-'' Clearly, at a given distance r*, the effect of the odd 
componc^nt in the Coulomb wave function is of increasing importance with a 
decreasing Bohr radius of the particle pair, i. e. for particles of greater masses 
or electric charges. At low energies, the sensitivity of the correlation to the 
odd component can be somewhat modified due to Coulomb interaction with 
the residual charge.^*^ 

It is clear that in the case of a dominant time asymmetry, v\{t)\ » K''!,)!, 
a straightforward way to determine the mean time delay (t) is to measure the 
correlation functions -R+(k*v > 0) and i?_(k*v < 0). Depending on the sign 
of (t), their ratio R+/R- should show a peak or a dip in the region of small 
k* and approach 1 at large values of k* . As the sign of the scalar product 
k*v is practically equal to that of the difference of particle velocities v\ — V2 
(this equality is always valid for particles of equal masses), the sensitivity of 
the correlation functions i?+ and i?_ to the sign of the difference of particle 

■'^This factor substantially deviates from unity only at k* < 2-ir/\a\ {e.g., at k* < 22 MeV/c 
for two protons). Note that for the distances r* > \a\ the confluent hypergeometrical function 
becomes important and compensates the deviation of the Coulomb factor from unity except 
for the classically forbidden region of k* < (|a|r*/2)~^/-^, narrowing with the increasing r*. 
Thus the FSl practically vanishes if at least one of the two particles comes from a long lived 
source {r],r]' , A, K'^, . . .). 



4 



emission times has a simple classical explanation. Clearly, the interaction 

between the particles in the case of an earlier emission of the faster particle 
will be weaker compared with the case of its later emission (the interaction 
time being longer in the latter case leading to a stronger correlation). This 
expectation is in accordance with Eqs. (1) and (4) at fc* ^ 0, (r*) ^ \a\ and 
(r*)|) < 1, when (the arrow indicates the limit v\{t) \ > |(rL)|): 

R^/R. « 1 + 2M ^ 1 _ 2^^^^^^-^^)^ (5) 
a a 

The sensitivity of the R+/R- correlation method to the mean relative 
time shifts (introduced ad hoc) was studiccP for various two-particle systems 
simulated in Pb+Pb collisions at SPS energy using the event generator VENUS 
5.14.^^ The scaling of the effect with the space-time asymmetry and with the 
inverse Bohr radius a, indicated by Eq. (5), was clearly illustrated for the 
K~^K~ system (a = —110 fm) and for the like- and unlike-sign nK, np and 
Kp systems (a = ±249, ±226 and ±84 fm respectively). It was concluded that 
for sufficiently relativistic pairs (71; > 0.5) the ratio can be sensitive to 

the shifts in the particle emission times of the order of a few fm/c. Motivated 
by this result the J?+/J?_ method was recently applied to the K^K~ system 
simulated in a two-phase thcrmodynamical evolution model and the sensitivity 
was demonstrated to the production of the transient strange quark matter state 
even if it decays on strong interaction time scales.^^ 

The method sensitivity was also studied for AGS and SPS energies using 
the transport code RQMD v2.3." Thus at SPS energy the central Pb + Pb 
collisions have been simulated and the 7r+i^+, 7r+p and K^p correlations have 
been studied.-'^^ To get rid of the effect of a fast longitudinal motion, the study 
was done in the longitudinally co-moving system (LCMS) in which the pair is 
emitted transverse to the reaction axis so that v = v_\_, = Ax and 

Ax* = 7_l(Ax - w^At), Ay* = Ay, Az* = Az. (6) 

The simulated correlation functions i?_ and their ratios are plotted in 
Fig. 1. We can see that for Tr+p and 7r+if+ systems these ratios arc less 
than unity at small values of q = k*, while for K'^p system the ratio i?+ 
practically coincides with unity. These results well agree with the mean values 
of At, Ax and Ax* presented in Table 1 {{Ay) « {Az) « 0). It can be 
seen from Eqs. (5) and (6) that the absence of the effect in the ratio 
for the if +p system is due to practically the complete compensation of the 
space and time asymmetries leading to Ax* « 0. For vr+p system the effect 
is determined mainly by the x- asymmetry. For n^K^ system both the x- 
and time-asymmetries contribute in the same direction, the latter contribution 
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being somewhat larger. The separation of the relative time delays from the 

spatial asymmetry is, in principle, possible (see Eq. (6)) by studying the ratio 
R+/R- in different intervals of the pair velocity. 

At AGS energy the Au + Au collisions have been simulated and the 7r+p 
correlations have been studied in the projectile fragmentation region where 
proton directed flow is most pronounced and where the proton and pion sources 
are expected to be shifted relative to each other both in the longitudinal and 
in the transverse (flow) directions in the reaction plane. It was demonstrated 
that a modification of the i?+ / R- method (± corresponding now to the signs 
of the respective components k*) is sufficiently sensitive to reveal these shifts.^^ 

At low energies, the particles in heavy ion collisions are emitted with the 
characteristic emission times of tens to hundreds fm/c so that the observ- 
able time shifts should be of the same order In fact the R+/R- method 
has already been successhilly applied to study proton-deuteron correlations in 
several heavy ion experiments at GANIL.^^'^^'^^ It was observed, in agreement 
with the coalescence model, that deuterons are on average emitted earlier than 
protons. 



4 Conclusion 

We have shown that unlike particle correlations, compared with those of iden- 
tical particles, contain a principally new piece of information on the relative 
space-time asymmetries in particle emission, thus allowing, in particular, a 
measurement of the mean relative delays in particle emission at time scales 
as small as 10~^^ s. To determine these asymmetries, the unlike particle cor- 
relation functions i?+ and i?_ have to be studied separately for positive and 
negative values of the projection of the relative momentum vector in pair c.m.s. 
on the pair velocity vector or, generally, - on any direction of interest. We have 
presented here the results of recent studies of these correlation functions for 
a number of two-particle systems simulated with various event generators. It 
was shown that the / i?_ ratio is sufficiently sensitive to the relative space- 
time asymmetries arising due to the formation of the quark-gluon plasma and 
strangeness distillation and even to those expected in the usual dynamical sce- 
narios at AGS and SPS energies. As to the detection of the unlike particles 
with close velocities (pi/mi « P2/'m2), there is no problem with the two-track 
resolution since these particles, having either different momenta or different 
charge-to mass ratios, have well separated trajectories in the detector mag- 
netic field. For the same reason, however, a large momentum acceptance of 
the detector is required. 



6 



Acknowledgments 

This work was supported by GA AV Czech RepubUc, Grant No. A1010601 
and by GA Czech RepubUc, Grant No. 202/98/1283. 

References 

1. G. Goldhaber et al., Phys. Rev. 120, 300 (1960). 

2. G.I. Kopylov, M.I. Podgorctsky, Yad. Fiz. 15, 392 (1972) {Sov. J. Nucl. 
Phys. 15, 219 (1972)); G.l. Kopylov, Phys. Lett. B 50, 472 (1974); 
M.I. Podgorctsky, Fiz. Elem. Chast. Atom. Yad. 20, 628 (1989) (Sov. 
.J. Part. Nucl. 20. 266 (1989)). 

3. R. Hanbury-Brown, R.Q. Twiss, Phil. Mag. 45, 663 (1954); Nature 178, 
1046 (1956). 

4. N.S. Amelin, R. Lednicky, Heavy Ion Physics 4, 241 (1996); SUBATECH 
95-08, Nantes 1995. 

5. S. Pratt, Phys. Lett. B 301, 159 (1993); Phys. Rev. C 50, 469 (1994). 

6. B. Erazmus ct al.. Internal Note ALICE 95-43. 

7. S.E. Koonin, Phys. Lett. B 70, 43 (1977); M. Gyulassy, S.K. Kauffmann, 
L.W. Wilson, Phys. Rev. C 20, 2267 (1979). 

8. R. Lednicky, V.L. Lyuboshitz, Yad. Fiz. 35, 1316 (1982) [Sov. J. Nucl. 
Phys. 35, 770 (1982)); Proc. CORINNE 90, Nantes, France, 1990 (ed. 
D. Ardouin, World Scientific, 1990) p. 42; JINR report P2-546-92 (1992); 
Hea,vy Ion Physics 3, 93 (1996). 

9. D.H. Boal, J.C. Shillcock, Phys. Rev. C 33, 549 (1986); D.H. Boal, C.-K. 
Gelbke, B.K. Jennings, Rev. Mod. Phys. 62, 553 (1990). 

10. R. Lednicky, V.L. Lyuboshitz, B. Erazmus, D. Nouais, Phys. Lett. B 
373, 30 (1996); Report SUBATECH 94-22, Nantes 1994. 

11. K. Werner, Phys. Rep. 232, 87 (1993); K. Werner and J. Aichelin, Phys. 
Rev. C 52, 1584 (1995). 

12. S. Soff et al., J. Phys. G12, 2095 (1997). 

13. H. Sorge et al., Phys. Lett. B 243, 7 (1990). 

14. R. Lednicky, S. Panitkin, Nu Xu, submitted to QM97, Tsukuba. 

15. S. Voloshin, R. Lednicky, S. Panitkin, Nu Xu, Phys. Rev. Lett. 79, 4766 
(1997). 

16. C. Ghisalberti et al., Nucl. Phys. A 583, 401 (1995). 

17. B. Erazmus et al., Proc. XXXIV Bormio Meeting (1996); D. Nouais, 
PhD Thesis (1996), Nantes (unpublished). 

18. D. Gourio, PhD Thesis (1996), Nantes (unpublished). 



7 



Table 1: Mean values of the relative space-time coordinates in LCMS (in fm) calculated 
from RQMD {v2.3) for 158 A-GeV Pb + Pb central collisions. 
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Figure 1: Unlike particle correlation functions i?+ and R— and their ratios simulated with 
RQMD for mid-rapidity particle pairs Tr+p, n'^K^, and X+p in Pb + Pb collisions at SPS 
energy. 
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